Giant hole and circular superflow in a fast rotating Bose-Einstein condensate 
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A fast rotating Bose-Einstein condensate confined in a quadratic-plus-quartic potential is found to 
dynamically generate a "giant vortex" that absorbs all phase singularities into a central low density 
hole, thereby sustaining a quasi-one-dimensional circular superflow at a supersonic speed. 
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Gaseous Bose-Einstein condensates (BECs) are an ex- 
tremely versatile testing ground for superfluidity particu- 
larly when an externally driven rotation exists. Through 
a number of experimental jl], |J and theoretical stud- 
ies H on rotating BECs, the argument now extends to 
the fast rotating regime || [?], ||, 0. For a harmonic 



trapping potential (l/2)muj 1 r , the centrifugal force pre- 
vents a BEC from rotating at frequency beyond the 
radial trapping frequency lo±. However, a recent exper- 
imental development jflj should enable one to create a 
confinement potential tighter than harmonic, thus open- 
ing a possible method to explore the nature of fast rotat- 
ing BECs. 

In this article, we extend our previous studies jl2j to 
the dynamics of vortex lattice formation of a BEC in a 
quadratic-plus-quartic potential with f2 greater than u)j_ . 
Here, the stable vortex configuration, realized at a value 
of Q much greater than u> ±_ , is not triangular but rather a 
circular array around a central low density hole [see Fig. 
1(e) below]. This hole is created where singly quantized 
vortices are close together but do not completely overlap. 
Surprisingly, as fl is further increased, all the vortices are 
absorbed into the central low density hole, around which 
persistent current circulates. This superflow is super- 
sonic, which offers an avenue of research in superfluidity, 
because it has been difficult to obtain a circular super- 
flow in superfluid helium |0|. For an atomic-gas BEC, a 
circular superflow can occur in an ideal situation and the 
superfluid speed is easily controlled by changing 0. This 
opens up possibilities of studying phase slippage (m) and 
the stability of a controlled superflow. 

We consider a two-dimensional system subject to ro- 
tation J7 = f2z by assuming translation symmetry along 
the z axis. In a frame rotating with frequency Q around 
the z axis, the dynamics of a condensate "wave function" 
^> is described by the Gross-Pitaevskii equation 



-V 2 + y trap -/i + C|<F| 2 -m z *, (1) 



where /u, is the chemical potential and the wave function 
is normalized to unity: J dxdy\^>\ 2 = 1. The units of 
energy, length, and time are given by the correspond- 
ing scales of the harmonic potential, that is, Huj±, b± = 
y / h/2muj±, and wj 1 , respectively, where m is the atomic 
mass and lu±_ the frequency of the harmonic potential. 



The strength of the mean- field interaction between atoms 
is characterized by C which is proportional to the s wave 
scattering length a s as C = 8irN2r>a s , where -/V 2 d is the 
particle number per unit length along the z axis. The 
trapping potential is Vt rap (r) = r 2 /4 + /cr 4 /16, where 
r 2 = x 2 + y 2 and k is the strength of the quartic term. 
The quartic part of Vt rap allows us to increase f2 above 
lu±. For a rotating condensate an effective potential in 
the radial direction is 
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(d). The lower 



FIG. 1: Time development of the angular momentum per 
atom l z (a) and the condensate density |\? 
(e). Times u ± t = 42 (b), 50 (c), and 78 
panel shows the density (e) and phase (f) profiles of the final 
stationary state of this dynamic system. The value of the 
phase varies continuously from (black) to 2ir (white). The 
discontinuity lines between black and white are the branch 
cuts in the complex plane. The termination of each branch 
cut is a vortex, or equivalently a phase defects. 
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This potential has a "Mexican hat" structure for > 1, 
which changes the dynamics of vortex lattice formation 
as shown below. 

Now, we investigate the dynamics of vortex lattice for- 
mation [fL5[| . The stationary solution of Eq. (0) for a 
nonrotating trap is used as the initial state. The rotat- 
ing drive is suddenly turned on by introducing a small 
anisotropy of the harmonic trap as (e x x 2 + e y y 2 )/4: pi, fL2[. 
Because the transition from a nonvortex state to a vor- 
tex state requires energy dissipation, a phenomenologi- 
cal damping constant 7 is put in Eq. (l|) by replacing 
id/dt with (i - j)d/dt @ [u|. Finally, the small trap 
anisotropy t XiV is turned off adiabatically to obtain the 
axisymmetric stationary state. 

Numerical simulations were done for several values of f2 
at fixed C = 250 and k = 1. Dynamics that are unique to 
the quadratic-plus-quartic potential occur when Q > 1. 
Figure shows the time development of the angular mo- 
mentum per atom l z = J dxdy^* L z ty , and the density 
and phase profiles for = 2.5. The ripples are excited 
on the surface [Fig. 0(b)], and some ripples penetrate 
into the condensate and develop into vortex cores [Fig. 
0(c)]. Some penetrating vortices move toward the center, 
merging together to make a hole at the center [Fig. 0(d)]. 
During the period between that shown in Fig. 0(c) and 
that in Fig. l](d) , additional vortices entered the conden- 
sate successively, which gradually increased the angular 
momentum. The density and phase profiles of the final 
steady state are shown in Figs. 0(e) and 0(f), respec- 
tively. A hole appears at the center, around which some 
vortices form a circular array. Furthermore, near the cen- 
ter of the condensate, some phase defects come very close 
to each other without overlapping [Fig. 0(f)]. 

As fl is further increased, all vortices generated from 
the condensate surface are absorbed into the central low 
density hole as shown in Fig. 0(a). Note that this central 
hole is composed of 24 singly quantized vortices packed 
together [Fig. g(b)]. This packing is possible for high ro- 
tation frequencies (O > 1) because the centrifugal force 
decreases the condensate density near the center; thus 
packing vortices costs little energy. The minimum of Eq. 
(pf) determines the radius of the ring condensate in Fig. 
g as R — -y/2(£l 2 — 1)1 k. Here we call a set of vortices 
such as that in Fig. @ a "giant vortex" to indicate that 
a number of phase defects are contained in a single hole 
[see Fig. 0(b)] . Such a density depression with a concen- 
tration of vorticity occurs in rotating superfluid helium 

Figures 0(f) and 0(b) show that the phase singularities 
do not completely overlap. The condensate inside the 
hole has a very low density and hence a long coherence 
length, so that the vortex cores have a nearly uniform 
density. Thus, only the locations of singularities matter 
and one can assume that the interaction energy between 
the singularity points is proportional to —pn 2 \nd, where 
p is the condensate density, k the quantum circulation, 
and d the distance between defects. The logarithmic di- 
vergence at d — > prevents the phase singularities from 
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FIG. 2: The condensate density (a) and phase field (b) of 
the vortex state with circulation 24/i/ra for = 3.2. The 
phase field region in (b) is the central low density region of 
(a), where 24 phase defects are marked by open circles. 



merging together. 

The appearance of such a density hole is relevant to re- 
cent studies on fast rotating BECs that used the Thomas- 
Fermi and Wigner-Seitz approximation [0, gj. Fischer 
and Baym evaluated the vortex core size self-consistently 
by minimizing the Gross-Pitaevskii energy functional and 
using a variational wave function they showed that 
vortex cores never overlap, in agreement with our finding. 
They also showed that the giant vortex state is energet- 
ically favorable at large il. Although the Thomas-Fermi 
approximation accurately describes stationary solutions 
such as that in Fig. 0(e), it fails when the width of the 
ring condensate becomes thin such as that in Fig. 0. 
Even if Q is not large enough, a pinning potential might 
be used to stabilize the giant vortex Jig] , 

For sufficiently large f2, the ring structure of Fig. 0(a) 
offers some insights into circular superflow, which has 
been studied for a torus pj. The superfluid velocity is 
given by the gradient of the phase of the wave function 
as v s = (h/m)V0; the circulation along a closed loop 
is quantized as <fv s ■ dl = nn with k = h/m being the 
quantum circulation. The quantum number n is equal 
to the number of branch cuts of the phase field along 
the loop or, equivalently, the number of phase defects 
in the low density hole. If there is sufficient dissipation, 
some vortices in the low density hole should spiral out of 
the condensate when Q is abruptly decreased from that 
of Fig. 0(a). The motion of these vortices across the 
circular superflow causes v s to decay. This is known as 
phase slippage which is a characteristic feature of 
superfluidity. 

Measurement of the angular momentum of the giant 
vortex allows the phase slippage to be observed well un- 
der control. For example, phase slippage was observed 
in helium superflow through microapertures |i~9|| , but it 
was difficult to control single phase slips and giant slips. 
Figure shows the time development of the angular mo- 
mentum and phase field resulting from a decrease in f2 
from = 3.5 (n = 32) to VI = 3.2. Here we assumed a 
small anisotropy in the trapping frequencies to break the 
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rotational symmetry. Some phase defects are released 
from the low density hole, jostling at the inner border 
of the ring condensate. We see that four vortices cross 
the superflow by phase slippage and then escape outside 
(only one escape is shown in Fig. ||) , thus decreasing the 
angular momentum. 




where 
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with A = J rdr\f\ 2 . The wave function is linearized as 
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resulting eigenvalue problem with the ansatz u(8) 



Uge i(n+q)e and v (g) 
relation 
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yields the dispersion 



(6) 




where q/R is the wave number of the excitation. The 
second term on the right-hand side of Eq. (Q) gives the 
sound velocity c = y/2CP. The superflow is stable when 
v s < c according to the Landau criterion. The values of 
the parameters and the form of f(r) of Fig. ||(a) yield 
c = 3.75 and the circular flow velocity v s = 2n/R ~ 12 
(in units of b±uj±). Because v s > c, the circular flow is 
supersonic and energetically unstable. 



FIG. 3: The top figure shows the time evolution of the angular 
momentum per atom £ z for a change of f2 from 3.5 to 3.2. The 
bottom figures show the phase profiles and line contours of the 
density for half of the x — y plane at cj±t = 102, 107, 112, and 
116 from left to right. Arrows show how one vortex escapes 
via phase slippage. 

When the ring radius R is sufficiently large, the circu- 
lar flow becomes quasi-one-dimensional for the following 
reason. The radial wave function of the stationary giant 
vortex state is obtained by substituting 5* = f(r)e m9 
into Eq. (|): 



(a) 



(b) 



d 2 1 8 n 2 

+ — + Krap + C\f\ 2 - HQ 

or z r or r z 



f=l*f. (3) 



The excitation of the radial component is easily obtained 
by writing / + Sf; the effective potential for 5f is then 
given by 



V cS (r) = Krap + -z + 2C\f{r)\ 2 -n- nfl. 



(4) 



Because n = flR 2 /2 in the limit of a rigid body ro- 
tation H 0], the radial zero-point energy around the 
minimum of the potential Vtrap + n 2 /r 2 becomes larger 
than the mean-field interaction strength C\f(r)\ 2 for 
large R. Such a high excitation gap freezes the mo- 
tion of the radial wave function while keeping the mo- 
tion free along the 9 direction. To investigate the sta- 
bility of this one-dimensional superflow, we substitute 
^{r,Q,t) = f(r)ip(6,t) into Eq. M to obtain 
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FIG. 4: Snapshots of the numerical simulation of a giant vor- 
tex state (Fig. ^ that is perturbed by an impurity potential. 
Parameter values for Vi mp are Vo = 100, a = 0.1, xo — 3.94 
and j/o = 0. The impurity potential is at rest in the laboratory 
frame, and thus rotates in the rotating frame. 

Impurities and walls can cause dissipation for a su- 
personic superflow. For example, Ref. |E| used the 
Gross-Pitaevskii equation for bulk condensates and or- 
dinary trapped condensates to show that the dissipa- 
tion of the superflow is caused by vortex-pair creation 
at an impurity object. This was shown experimentally 
pl[ . However, this does not apply to the present case 
of quasi-onc-dimensional superflow because of the tight 
radial confinement p2[ . To show the influence of impuri- 
ties, we introduce an impurity potential of Gaussian form 
as V; 



imp 



V r exp[-{(a;-xo) 2 + (y-yo) 2 }/o' 2 ] (at rest 



the laboratory frame), where {xq, j/o} is the impurity po- 
sition. Figure ^ shows an example of the destabilization 
process of the giant vortex state by Vimp, for which Vq 
is chosen to be larger than T4rap(-R) an d & smaller than 
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the width of the ring condensate. Because v s exceeds the 
Landau critical velocity c, the impurity excites the lon- 
gitudinal density waves along the ring condensate, thus 
producing a drag force. Due to the presence of dissipation 
and the rotating drive, the system reaches a steady state 
as shown in Fig. (|(b). Therefore, the quadratic-plus- 
quartic system offers a scheme for studying the stability 
of a quasi-one-dimensional superflow. 

In conclusion, a fast rotating BEC in a quadratic-plus- 
quartic potential can generate a giant vortex and cir- 
cular superflow around it by absorbing singly quantized 
vortices into a central low density hole. The circular 



superflow is quasi-one-dimensional and goes at a super- 
sonic speed. We thus hope that this system will encour- 
age experiments to study low dimensional superfluidity 
in atomic-gas BECs. 
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